Multidimensional gravity in non-relativistic limit 
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It is found the exact solution of the Poisson equation for the multidimensional space with topology 
M'i+d = K 3 x T d . This solution describes smooth transition from the newtonian behavior l/r3 for 
distances bigger than periods of tori (the extra dimension sizes) to multidimensional behavior l/Vgijj 
in opposite limit. In the case of one extra dimension d — 1, the gravitational potential is expressed 
via compact and elegant formula. It is shown that the corrections to the gravitational constant in 
the Cavendish-type experiment can be within the measurement accuracy of Newton's gravitational 
constant Gn- It is proposed models where the test masses are smeared over some (or all) extra 
dimensions. In 10-dimensional spacetime with 3 smeared extra dimensions, it is shown that the size 
of 3 rest extra dimensions can be enlarged up to submillimeter for the case of ITeV fundamental 
Planck scale M P i( W y In the models where all extra dimensions are smeared, the gravitational 
potential exactly coincides with the newtonian one. Nevertheless, the hierarchy problem can be 
solved in these models. 
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Introduction There are two well-known problems 
which are related to each other. They are the discrepan- 
cies in gravitational constant experimental data and the 
hierarchy problem. Discrepancies (see e.g. Figure 2 in 
the "CODATA Recommended Values of the Fundamen- 
tal Constants: 2006") are usually explained by extreme 
weakness of gravity. It is very difficult to measure the 
Newton's gravitational constant Gn- Certainly, for this 
reason geometry of an experimental setup can effect on 
data. However, it may well be that, the discrepancies 
can also be explained (at least partly) by underlying fun- 
damental theory. Formulas for an effective gravitational 
constant following from such theory can be sensitive to 
the geometry of experiments. For example, if correction 
to the Newton's gravitational potential has the form of 
Yukawa potential, then the force due to this potential 
at a given minimum separation per unit test-body mass 
is least for two spheres and greatest for two planes (see 
e.g.[3|). The hierarchy problem - the huge gap between 
the electroweak scale Mew ~ 10 3 GeV and the Planck 
scale Mp;(4) = 1.2 x 10 19 GeV - can be also reformulated 
in the following manner: why is gravity so weak? The 
smallness of Gat is the result of relation Gn = Mr 
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The natural explanation was proposed in 2, 3]: the grav 



ity is strong: Gt> = M 



(2+d) 
Pl{T>) 



M 



-(2+d) 
EW 



and it happens 



in (2? = 4 + d)-dimensional spacetime. It becomes weak 
when gravity is "smeared" over large extra dimensions: 
Gn ~ G-p/Vd where Vd is a volume of internal space. 
To shed light on both of these problems from new stand- 
point we intend to investigate multidimensional gravity 
in non-relativistic limit. 

Multidimensional gravitational potentials It is of in- 
terest to generalize the well-known Newton's gravita- 
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tional potential ^p{r^} = — Gjvm/r3 {r§ — \r^\ is mag- 
nitude of a radius vector in three-dimensional space) to 
multidimensional case. Clearly, the result depends on 
topology of investigated models. We consider models 
where (23 — 3+d)-dimensional spatial part of factorizable 
geometry is defined on a product manifold Mr) — M. 3 xT d . 
K 3 describes three-dimensional flat external (our) space 
and T d is a torus which corresponds to d-dimensional 
internal space with volume Vd- Let b ~ Vj be a char- 
acteristic size of extra dimensions. Then, Gauss's flux 
theorem leads to the following asymptotes for gravita- 
tional potential (see e.g. [H): <p ~ l/r 3 for r 3 >> b and 
tp ~ 1/ r^t d for r% + d << b where r% + d is magnitude of a 
radius vector in (3 + (i)-dimensional space. 

To get the exact expression for 23-dimensional gravita- 
tional potential, we start with the Poisson equation: 



& D ip D = S d Gt>Pd{*d) , 



(1) 



where Sd = 2tt d ^ 2 /T(D/2) is a total solid angle (square 
of (D — l)-dimensional sphere of a unit radius), Gp is a 
gravitational constant in {T> = D + l)-dimensional space- 
time and pd{vd) — mS(xi)5(x2)--S(xr))- In the case of 
topology R D , Eq. ([T]) has the following solution: 
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23 > 3. 
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This is the unique solution of Eq. (p} which satisfies 
the boundary condition: lim <pd{yd) = 0. Gravita- 

rn— >+oo 

tional constant Gr> in Jlj) is normalized in such a way 
that the strength of gravitational field (acceleration of a 
test body) takes the form: —dtpo/dro = -Gpm/r^ 1 . 

If topology of space is M 3 x T d 1 then it is natural to 
impose periodic boundary conditions in the directions 
of the extra dimensions: <^r>(r3, £i, £21 • • • , £i, ■ ■ ■ , £d) = 
<PD(r3,€i,&, ■ ■ ■ ,€i + a n ■ ■ ■ i = l,...,d, where ai 

denotes a period in the direction of the extra dimension 
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Then, Poisson equation has solution (cf. also with 
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To get this result we, first, use the formula <5(£i) = 
~ 5Zfe^_oo cos (iTT&J an( i' second, put the following re- 
lation between gravitational constants in four- and T>- 
dimensional spacetimes: 



Sd 
S 3 
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The letter relation provides correct limit when all eij — > 0. 
In this limit zero modes fcj = give the main contribution 
and we obtain ^n(r 3 , £1, ^) — > ~GNm/r 3 . Eq. (gj 
was widely used in the concept of large extra dimensions 
which gives possibility to solve the hierarchy problem 0, 
3]. It is also convenient to rewrite ^ via fundamental 
Planck scales: 



— • M 2 - M 2+d 

g m Pl(i) — IVI Pl(V) 
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where M 



G 



-l/(2+d) 



Pi(4) 



G 



-1/2 



/V 



1.2 x 10 iJ GeV and M PiCD) 



p 



are fundamental Planck scales in four and T> 
spacetime dimensions, respectively. 

In opposite limit when all — ► +00 the sums in Eq. 
([3]) can be replaced by integrals. Using the standard 
integrals (e.g. from |5j) and relation we can easily 
show that, for example, in particular cases d — 1, 2 we get 
desire result: <££>(r 3) £x, ■■■■>€d)-+ -G v m/[(D-2) r\~l% 
One extra dimension In the case of one extra dimen- 
sion d = 1 we can perform summation of series in Eq. 
((3|). To do it, we can apply the Abel-Plana formula or 
simply use the tables of series [5]. As a result, we arrive 
at compact and nice expression: 



¥>4(>3,£) = 
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(6) 



where r% E [0, +00) and £ E [0, a]. It is not difficult 
to verify that this formula has correct asymptotes when 
r.3 >> a and r± << a. Fig. 1 demonstrates the shape 
of this potential. Dimensionlcss variables rji = r^/a E 
[0, +00) and 772 = £/a E [0, 1]. With respect to variable 
772, this potential has two minima at 772 = 0, 1 and one 
maximum at 772 = 1/2. We continue the graph to nega- 
tive values of 772 E [— 1, 1] to show in more detail the form 
of minimum at 772 = 0. The potential © is finite for any 



value of r3 if £ ^ 0, a and goes to —00 as —\jr\ if simul- 
taneously 7-3 — ► and £ — + 0, a (see Fig. 2). We would 
like to mention that in particular case £ = formula ([6|) 
was also found in Q. 




FIG. 1: Graph of function ^(771,772) = V?4(r3, £) / (GNm/a) 
— sinh(27r77i)/[77i(cosh(27rr;i) — cos(27T772))]. 
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FIG. 2: Section £ = of potential $6$. Solid line is ^(771, 0) = 
— sinh(27r77i)/[77i (cosh(27T7;i) — 1)] which goes to — 1/771 (dot- 



ted line) for ni 
771 0. 



+00 and to — 1/(tt77i) (dashed line) for 



Having at hand formulas ([3]) and ([5]), we can apply it 
for calculation of some elementary physical problems and 
compare obtained results with known newtonian expres- 
sions. Some of these calculations can be found in our 
preprint Q. For a working approximation, it is usually 
sufficient to summarize in ([3|) up to the first Kaluza-Klein 
modes |fcj| = 1 (i = 1, . . . , d). Then, the terms with the 
biggest periods aj give the main contributions. If all test 
bodies are on the same brane (£, = 0) we obtain: 



¥>D(r3,6 = 0, = 0) = (p D (r3) 
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T3 



1 + a exp 
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(7) 



where a = 2s (1 < s < d) , A = a/(2n) and s is a 
number of extra dimensions with periods of tori eij which 
are equal (or approximately equal) to a — maxa^. If 
ai = CI2 = • • • = a,d = a>, then s = d. Thus, the correction 
to Newton's potential has the form of Yukawa potential. 
It is now customary to interpret tests of gravitational 
inverse-square law (ISL) as setting limits on additional 
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Yukawa contribution. The overall diagram of the experi- 
mental constraints can be found in [lj (see Figure 5) and 
we shall use these data for limitation a for given a. 

In this approximation, the gravitational force between 
two spheres with masses mi, ma, radiuses R%, R 2 and dis- 
tance r 3 between the centers of the spheres reads: 



F 



G N ( e ff)m 1 m 2 



(8) 



where 



2irr 3 
x exp 



2tt 



(r 3 ~R 1 -R 2 ) 



2nR 2/ 
G N (1 + 6 G ).(9) 



Here, we used conditions: r 3 > Ri + R 2 and Ri,R 2 
a/2n. 

Smeared extra dimensions In what follows, we con- 
sider asymmetrical extra dimension model (cf. Q) with 
topology 



M D = K 3 x T d ~ p xf, p < d , 



(10) 



where we suppose that {d — p) tori have the same "large" 
period a and p tori have " small" equal periods b. In this 
case, the fundamental Planck scale relation ([5]) reads 
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Additionally, we assume that test bodies are uniformly 
smeared/ spreaded over small extra dimensions. Thus, 
test bodies have a finite thickness in small extra dimen- 
sions (thick brane approximation). For short, we shall 
call such small extra dimensions as "smeared" extra di- 
mensions. If p = d then all extra dimensions are smeared. 

It is not difficult to show that the gravitational po- 
tential does not feel smeared extra dimensions. We can 
prove this statement by three different methods. First, 
we can directly solve D-dimensional Poisson equation 
([T|) with the periodic boundary conditions for the ex- 
tra dimensions £ p+ i,...,£d and the mass density p — 
( TO /nLi a.i) 8(r 3 )6(€ p+1 )...6(€d)- Second, we can aver- 
age solutions ([3]) and ([6]) over dimensions £1,.. -i^p an d 
take into account that J Q a cos(27rfc£/a)<i£ = 0. In par- 
ticular case of one extra dimensional, we can also show 
that 
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(12) 



Finally, it is clear that in the case of test masses smeared 
over extra dimensions, the gravitational field vector 
~Ed = —^d^Pd does not have components with respect 



to extra dimensions: E_d = Eoti r3 . Thus, applying the 
Gauss's flux theorem to the Poisson equation, we obtain: 
Eo{r 3 ) = -G N m/r% -> <££>(r 3 ) = -G N m/r 3 . There- 
fore, all these 3 approaches show that in the case of p 
smeared extra dimensions the wave numbers ki , . . . , k p 
disappear from Eq. Q and we should perform summa- 
tion only with respect to fc p +i, . . . , kd- 

Effective gravitational constant Coming back to the 
effective gravitational constant ([9]) in the case of topol- 
ogy (fT0|) with p smeared extra dimensions, s in Eq. (|9|) 
is replaced by (d — p). Now, we want to evaluate the cor- 
rections Sq to the Newton's gravitational constant and 
to estimate their possible influence on the experimen- 
tal data. As it follows from Figure 2 in the CODATA 
2006, the most precise values of Gn were obtained in 
the University Washington and the University Zurich 
experiments [| They are Gjv/10 _u m 3 ks _1 s- 2 = 

6.674215 ± 0.000092, and 6.674252 ± 0.000124, respec- 
tively. Let us consider two particular examples: (T> = 5)- 
dimensional model with d = l,p = — » a = 2 and 
(T> = 10)-dimensional model with d = 6,p — 3 — » a = 6. 
For these values of a, Figure 5 in [l[ gives the upper limits 
for A = a I (2ir) correspondingly A«2x 10~ 2 cm and A w 
1.3 x 10 _2 cm. To calculate Sg, we take parameters of the 
Moscow experiment [ll[: R\ ~ 0.087cm for a platinum 
ball with the mass mi = 59.25 x 10 _3 g, R 2 ~ 0.206cm 
for a tungsten ball with the mass m 2 = 706 x 10~ 3 g 
and r 3 — 0.3773cm. For both of these models we obtain 
S G w 0.0006247 and S G » 0.0000532, respectively. Both 
of these values are very close to the measurement accu- 
racy of Gn in @, 03 ■ So, if the same accuracy can be 
achieved in the Moscow-type experiments, then, chang- 
ing values of R\ y2 and r 3 , we can reveal extra dimensions 
or obtain experimental limitations on considered models. 

Model: V — 10 with d = 6,p = 3 Let us consider in 
more detail (T> — 10)-dimensional model with 3 smeared 
dimensions. Here, we have very symmetrical with respect 
to a number of spacial dimensions structure: 3 our exter- 
nal dimensions, 3 large extra dimensions with periods a 
and three small smeared extra dimensions with periods 
b. For b we put a limitation: 6 < b max — 10~ 17 cm which 
is usually taken for thick brane approximation. As we 
mentioned above, in the case of a = 6, for a we should 
take a limitation a < a max = 8.2 x 10~ 2 cm. To solve the 
hierarchy problem, the multidimensional Planck scale is 
usually considered from ITeV up to approximately 130 
TeV (see e.g. @, H3|)- To make some estimates, we take 



M, 



= 1TcV< M 



Pino) 



< 



M max = 50TeV. Thus, as 
it follows from Eq. (TPI) . the allowed values of a and b 
should satisfy inequalities: 



c M 2 

09 lu Pi(i) 

*~^0 1 11 v 



< a 3 b 3 < 



Sg M Pl(i) 
S 3 AfjL„ 



(13) 



Counting all limitations, we find allowed region for a and 
6 (shadow area in Fig. ([3])). In this trapezium, the up- 
per horizontal and right vertical lines are the decimal 
logarithms of a max and b max , respectively. The right 
and left inclined lines correspond to Mp;(i ) = ITeV 



4 



and Mpz(io) = 50TeV, respectively. To illustrate this 
picture, we consider two points A and B on the line 
Mpi(io) = ITeV. Here, we have a = 0.82 x lO^cm, 
b = 10 _21 - 5 cm for A and a = 10~ 4 cm, b = 10" 18 (3 cm 
for B. These values of large extra dimensions a are much 
bigger than in the standard approach a ~ 10 < - 32 / 6 - ) ~ :L7 cm 
S3 10- n - 7 cm Hi. 



lg a 




FIG. 3: Allowed region (shadow area) for periods of large 
(a) and smeared (b) dimensions in the model T> = 10 with 
d — 6,p — 3. 

Model: T>-arbitrary and d = p In this model, the test 
masses are smeared over all extra dimensions. There- 
fore, in non-relativistic limit, there is no deviation from 
the Newton's law at all. It worth of noting, that this 
result does not depend on the size of smeared extra di- 
mensions. The ISL experiments will not show any devi- 
ation from the Newton's law without regard to the size 
b (see also Eq. ([9]) where s = d — p = 0). Similar rea- 
soning are also applicable to Coulomb's law. It is neces- 
sary to suggest other experiments which can reveal the 
multidimensionality of our spacetime. Nevertheless, we 
can solve the hierarchy problem in this model because 
Eq. (where d = p) still works. For example, in 

the case of bosonic string dimension T> = 26 we find 
Mp;( 26 ) ~ 31TeV for b = lCr 17 cm. In the case V = 10 
we get M Pl{w) ks 30TeV for b = 5.59 x 10- 14 cm. 

Conclusions We have considered generalization of the 
Newton's potential to the case of extra dimensions where 



multidimensional space has topology Mo = M 3 x T d . 
We obtained the exact solution which describes smooth 
transition from the newtonian behavior 1 /r% for distances 
bigger than periods of tori (the extra dimension sizes) to 
multidimensional behavior l/r D ~ 2 in opposite limit. In 
the case of one extra dimension, the gravitational poten- 
tial is expressed via compact and elegant formula (jHJ). 

As an Yukawa potential approximation, it was shown 
that the corrections to the gravitational constant in the 
Cavendish-type experiment can be within the measure- 
ment accuracy of Gjv- It may reveal the extra dimen- 
sions or provide experimental limitations on parameters 
of multidimensional models. 

Then, we proposed models where test masses can be 
smeared over extra dimensions. In this case, the gravita- 
tional potential does not feel smeared dimensions. The 
number of smeared dimensions can be equal or less than 
the total number of the extra dimensions. Such approach 
opens new remarkable possibilities. 

For example in the case V = 10 with 3 large and 3 
smeared extra dimensions and Mp;( 10 ) = ITeV, the large 
extra dimensions can be as big as the upper limit es- 
tablished by the ISL experiments for a = 6, i.e. a « 
0.82 x 10 _1 cm. This value of a is in many orders of mag- 
nitude bigger than the rough estimate a w 10~ u cm 
obtained from the fundamental Planck scale relation of 
the form of Eq. ©. 

The limiting case where all extra dimensions are 
smeared is another interesting example. Here, there is 
no deviation from the Newton's law at all. Nevertheless, 
the hierarchy problem can be solved in this model. 
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